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Abstract 

The quasiparticle excitation spectrum of isolated vortices in clean layered <i-wave superconduc- 
tors is calculated. A large peak in the density of states in the "pancake" vortex core is found, 
in an agreement with the recent experimental data for high-temperature superconductors. 

1 Introduction 

The vortex core in classical type-II superconductors can be treated as the normal metal area 
with radius of the order of the coherence length £(0) ~ lOnm. |], § The spectrum of the bound 
states near the Fermi surface, formed by the constructive interference between the incident 
and the Andreev reflected quasiparticles, @ is quasicontinuous (gapless). High-temperature 
superconductors (HTS) are layered, having a cylindrical Fermi surface, the superconductivity is 
of the strong and c?-wave coupling type, and the vortex core radius is much smaller, £(0) ~ lnm. 

For a two-dimensional (2D) vortex and s-pairing, Rainer et al. have shown, using the Andreev 
quasiclassical theory in the analytical, and Eilenberger's in the numerical part of their study, 
that bound states exist in the vortex core.[|J Similar conclusions have also been obtained by 
Maki and coworkers for ci-pairing, within the Bogoliubov-de Gennes approach. 

In this paper, the Eilenberger quasiclassical equations || , in the case of both s and (i-pairing, 
are solved analytically, within the model considering the spatial variation of the order parameter 
in the vortex core as for a normal metal cylinder of radius r c ~ £. A crucial difference is found 
in the quasiparticle spectra below the bulk energy gap, between the classical superconductors 
with spherical Fermi surface, s-pairing, large £, and HTS with cylindrical Fermi surface, d- 
pairing, small £. Our results confirm that the quasiparticle density of states (DOS) has one 
large maximum, recently observed in YBCO by scanning tunneling microscopy. 



2 Model of the Vortex Core 

An efficient method for calculating local spectral properties is the quasiclassical theory of su- 
perconductivity, which gives the Eilenberger equations 
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Here g = g^(r,v,cu n ) and / = f^(r,v,u n ) represent the normal and the anomalous Green 
function respectively, A = A(r, v) is the gap function, uj n = irk B T{2n + 1), n = 0, ±1,±2..., 
are the Matsubara's frequencies, and v is the Fermi velocity vector. The function f< is defined 
by /ji(r, v, u n ) = /*|(r, —v,u n ). A and / are connected by the self-consistency equation. 

For a homogeneous and isotropic superconductor, solutions of the Eilenberger equations are 
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In the cylindrical coordinates (r,(p,z), with the origin situated on the vortex axis, the vortex 
magnetic field is h = he z . At distances r ~ ( <C A from the vortex axis, h is approximately 
constant, and the gauge can be chosen in the form A = (h x r) /2. Taking the same vortex 
gap function as for a normal metal cylinder embedded in a superconductor 
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Here, r c is the vortex core radius, and 9 is the polar angle in k-space. For c/-pairing A(6) = 
A cos 2$, and for s-pairing A(9) = A . The gauge in Eq. (f|) is due to the flux quantization. 

For a pancake vortex in (r, (p) plane, denoting the coordinate along v by s, and along h x v by 
p (Fig.l.), in the gauge with real gap, Eqs. (|I|) and (0) can be rewritten in the form 
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where r 2 = p 2 + s 2 and l 2 H = he/ 'eh. 

For a normal metal cylinder and zero magnetic field, Eqs. ©-(0) with p = 0, the solution is of 
the form 

/ = E/*(p) eW . 9 = T,9i(p)e KiS - (8) 



For r < r c , with k = 2u n /v, 
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For r > r c and k = 2e n /hv, 
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Using the continuity condition for / and g at ±sq = ±y ?" 2 — p 2 , for r < r c the normal Green 
function is 
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For a vortex, approximating p/r 2 by p/r 2 , the solution of Eq. ©-(0) can be obtained from Eq. 
(|13|), by changing u n 
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and = 2u' n /v. In this case, the magnetic flux quantization leads to [l], ||] 
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Since for an isolated vortex Ih 3> r c , the direct influence of the field can be neglected, and the 
only relevant contribution is due to the screening supercurrent flow, ipv/2r 2 term in Eq. ([15]). 



3 Bound States 

Performing an analytical continuation of G by hu n — > —iE+7], E being the quasiparticle energy 
with respect to the Fermi level, the retarded propagator g R (E,p,8) is obtained. In terms of 



reduced variables E/A -> E, ^JA 2 (9) - E 2 /A Q ->■ e, y E 2 — A 2 (8)/A -> e, p/r c -> p, 
2so/vr^ — >■ so, £o = hv/7iA being the BCS coherence length, angle resolved partial DOS 
(PDOS) is obtained from N(E,p,6) = Reg R (E,p,6). 

For the normal metal cylinder 
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where 5 is the Dirac function, is the step-function, and iV(0) = m/27rh 2 is the normal metal 
density of states at the Fermi surface for one spin orientation. For s-wave pairing, PDOS does 
not depend on 9, while for d-wave pairing, averaging over the cylindrical Fermi surface leads to 

N(E,p)/N(0) = — \ 2 " N(E,p,9)/N(0)d9 = 



2vr Jo 



E \ \E\e 2 de 

+ 



7T J^m^{o,Ei-i} V£ 2 - eVl — E 2 + e 2 (e 2 cos 2 (Es ) + E 2 sm 2 (E 



so 



+ (^tan(E g0 ) + |^tan(^ )|) e ( _ ^ 

/cos 2 (Es ) - -E 2 V 7 



Finally, after spatial averaging over the cylinder area 7rr 2 , DOS is 

N(E) = - f 1 N(E,p)Jl-p 2 dp. (19) 
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For the vortex, in Eqs. fllTD and flU), E —> E + Eq signE, £ = ft|pi|v/2A r 2 , Eq. (|TJ), with p; 
from Eq. (|T6|) , and J2 Pi instead of integration in Eq. (|19"D . Here, signs of p and .E are connected, 



because the magnetic field causes a difference in propagation of particles and holes. 



For small radius vortices in HTS, ^omv/H ~ 1 (~ 10 in classical superconductors), only one 
trajectory through the vortex core is allowed, with p = p , Eq. (|i~6l). Taking for YBCO r c = £ , 
p = 1/3, A /E F = 0.424, only one peak in DOS in the vortex core around E/A Q m 0.3 
is obtained (Fig. 2). For comparison, DOS of normal metal cylinder embedded in the same 
superconductor and with the same radius r c = £ , but in the zero magnetic field, is shown. In 
this case, a large energy gap is found in DOS, due to formation of lowest bound state at high 
energy, of the order of A . This is not the case in classical superconductors, where A /Ef <C 1. 



Fig.l. Trajectory passing at dis- 
tance p from the vortex center. 

Fig. 2. Quasiparticle DOS in the vortex core (solid 
curve), and in the normal metal cylinder (dashed curve), 
for clean layered <i-wave superconductor. r c /£ = 1, 
A /E F = 0.424. 

In conclusion, cylindrical Fermi surface, <i-wave pairing and small £o, large A /Ep ~ 0.1, make 
DOS of a normal metal cylinder embedded in HTS and a pancake vortex different from DOS of 
a normal cylinder and a vortex in classical superconductors. Since the Andreev bound states 
can transport charge currents, unlike the bound states in a potential well, supercurrents can 
flow through the vortex without losses, strongly influencing its dynamics. This could be very 
important for transport properties of HTS, especially for understanding the unusual magnetic- 
field dependence of the electrothermal conductivity, which was observed experimentally and 
awaits explanation. || 
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